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The contact value of the charge profile for nonsymmetrical electrolytes is presented as the sum 
of three contributions. One of them is the normal component of the Maxwell electrostatic stress 
tensor. The second one is the surface electrostatic property defined as the integral of the product 
of the gradient of the electrical potential and the density distribution function of colons. The third 
term is the bulk contribution defined by the sum for anions and for cations of the product of their 
charge and their partial pressure. For noncharged surfaces only the last two are present and have 
the same sign in the case of size asymmetry. In the case of charge asymmetry the contact value of 
the charge profile is the result of the competitions of bulk and surface terms in which the bulk term 
is dominant. Using both the contact theorems for the density and the charge profiles, the exact 
expressions for the contact values of the profiles of colons and counterions are obtained and some 
related properties are discussed. 

The properties of the electrical double layer are described by the density and charge profiles p{z) and q{z). They 
are connected with the density distribution functions Pa{z) for ions of type a having the charge by the relations 

p{z) = (1) 
q{z) = ^eaPa{z) (2) 

a 

where z is the normal distance from the wall. 

For the contact value (CV) of the density and charge profiles, for the primitive model of electrolyte solutions, exact 
relations have established known as the contact theorems (CTs). The CT for the density profile states that the CV 
of the number density profile p^* is given by the sum of the bulk electrolyte pressure P and the electrostatic Maxwell 
stress tensor contribution [ij, @| 

/*^p(d/2)=/3P + /3— (3) 

where d is the ion diameter, assumed to be the same for all species, /3 = l/(fcBT) with fc^ the Boltzmann constant 
and T the absolute temperature, e is the dielectric constant of the solvent, eE/An — qs is the surface charge density 
per unit area on the wall. The CT for the charge profile was formulated recently by Holovko et al. fSj. In contrast to 
the theorem for the density profile, the CV of the charge profile q'^* has a nonlocal character and in general can be 
presented in the form 



q-t^q(d/2)=p r p{z)^^dz + py2e^P^ (4) 
Jd/2 dz ^ 

where p[z) = ej^paiz), ^piz) is the electrical potential 

An f°° 

ip{z) = — — J q{zi){zi - z)dzi (5) 
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and Pa is the bulk partial pressure for ions of type a. 

Both CTs ([3]) and ([U can be obtained by direct integration of the system of Born- Green- Yvon (BGY) equations for 
the singlet distribution functions pa (z) [1, Q . This procedure can be easily repeated for an electrolyte with different 
ion diameters da and leads to changing only p^' and q^* in relations ^ and ^ by the corresponding CVs given by 

p^* = (6) 
q^' = ^e„p„K/2) (7) 

a 

and the bulk partial pressure Pa according to the virial expression for the bulk pressure of a ion hard sphere system 
is given by 



pCoul 



(8) 



where da^y — (da + d~f)/2, ga-yidaj) is the CV of the pair distribution function ga^yir) and 

-1 />oo 

E^"""' = 7^P<. E 47rp,^ / 9a,{r)rdT (9) 

is the coulombic part of the partial internal energy. 

For the symmetrical electrolyte p{z) = e^p(z) (e being the elementary charge), it was shown [6] that the CV of the 
charge profile can be presented as the sum of the normal component of the Maxwell electrostatic stress tensor and a 
new electrostatic property defined by the integral of the product of the gradient of the electrical potential and the 
singlet distribution function of coions (ions having the identical sign as the surface charge). In this note we generalize 
this result for nonsymmetrical electrolytes and discuss some specific properties connected with the charge and size 
asymmetry of cations and anions. 

First, we consider a nonsymmetrical electrolyte with a single type of anion and of cation, p±, e± and d± are 
respectively the densities, charges and diameters of the ions and p{z) can be written as 

p{z) = e\p+{z) + etp-{z) (10) 
= e=pg(z) - e±p±(z)(e^ - e±) (11) 

where the upper/lower sign are relative to a positively/negatively charged surface. As a result the CT for the charge 
profile for nonsymmetrical electrolyte can be presented in the following form 

where 

. r dz'-^p^iz) (13) 

Jd±/2 oz 

p±{z) is the density distribution function for coions near a positively/negatively charged surface. The result obtained 
is the generalization of our previous result [1] for the CV of the charge profile of the symmetrical electrolyte to the 
nonsymmetrical case. Similar as for the symmetrical electrolyte expression (|12p has a simple interpretation. The 
difference between forces from counterions and coions required to obtain g'^* is connected with three contributions: 
the electrostatic Maxwell stress tensor contribution, the electrostatic contribution from coions described by the elec- 
trostatic property ^+ or JTL depending on the sign of the surface charge and a special combination of the partial 
pressures for coions and counterions. 

From the CTs for the density and charge profiles presented in their form ^ and (H^), it is possible to obtain the 
following exact relations for the CVs of the density profiles for coions and counterions respectively 

pf = p+{d+/2) = pP++(3e+J+, 

p-1 = p-(d_/2) = /3P_ + /3— - pe+J+ (14) 
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for a positively charged surface or 

= p+{d+/2) - [3P+ + - /?e_ J„, 

p!* EE p_(d_/2) = /3P_+/3e^J_ (15) 

for a negatively charged surface. Since the CV of the density profile of colons cannot be negative, we also have the 
following conditions 

P+ + e+J+ > (16) 

for a positively charged surface or 

P-+e_JL>0 (17) 

for a negatively charged surface. 

These inequalities reduce to equalities in the limit E —> oo and as a consequence we have an exact cancellation 
of e+J+ (or e-J-) by the partial pressure P+ (or P-). Similar as for the symmetrical electrolyte in the limit of 
small -E, Jj^ (or J) has a linear dependence on E and we can conjecture that Jj^ (or J J) is a monotonic function of E. 



However, for a noncharged surface {E = 0) and nonsymmetrical electrolytes, the charge profile q(z) is not zero in 
contrast to the symmetrical case. The charge and/or size asymmetry of cations and anions induces a double layer 
with an associated surface potential 13,181]. According to the CT eq. ([T^ . the CV of the charge profile for a noncharged 
surface is defined by two terms 

g^* = -/3e+(e_-e+)J++/3(e+P+ + e_P_) (18) 
= -/3e_(e+-e_)J_+/?(e+P+ + e_P_). (19) 

Both forms for the CV of the charge are equivalent and due to this for a noncharged surface we have the following 
symmetrical relation 

e+JV + e^J^ =0. (20) 

We may now distinguish the electrical double layer due to the charge and size asymmetry of ions for they are 
different in nature. 

The influence of size asymmetry has a short-ranged steric origin since the smallest ion can approach the surface more 
closely than the larger ions. If we consider for simplicity a charge symmetrical electrolyte (e+ — — e_ — e) according 
to relation 

(7=* = 2/3e2j++/3e(P+-P_). (21) 

Due to electroneutrality, the first layer constituted of the smaller ions has to be compensated by the following layer 
where we find the larger ions. Then intuitively, the CV of the profile for the larger ions which are further from 
the wall has to be larger in order to compensate the profile for the smaller ions which starts before as they can 
come closer to the wall. For example, for the case d+ < c?_ we will have p^j^ < pi* and g"^* will be negative. This 
intuitive result is indeed in agreement with the CT for the charge given in eq. (PT|) . For this case [(hli[z) j dz) < and 
according to definition eq. J'+ < thus the first term of eq. (PT|) is negative. The second term at least for small 
ion concentration is also negative since the dominant contribution to the partial pressure comes from the Coulombic 
part of the partial internal energy eq. ([9]) which is more negative for smaller ions. This means that in this case g"^* 
will be negative in agreement with the intuitive statement. For the opposite case d+ > d^, p'^ > pf* and g"^* will be 
positive. This result also corresponds to the CT (pij) as in this case {dip{z)/dz) > 0, J+ > and in parallel P_|- > P_ 
and g'^* will be positive. 

The infiuence of the charge asymmetry on is different. It is associated with the electrostastics and results from 
the long-range behaviour of the profiles of the cations and anions leading to an exclusion effect by the neutral surface 
[l, 8]. From \^], whe n the conditions are such that electrostatic interactions become important in comparison to 
steric effects, we know for symmetrical electrolytes that the density profile decreases when we approach the wall 
producing a negative adsorption of ions. This effect is more pronounced the higher the charge of the ions. Now 
this phenomenon can be generalized to ions with different charges where to simplify we consider size symmetrical 
electrolytes {d+ — d^ ^ d). And as a consequence, the CV of the charge for the ion carrying the higher charge 
will be smaller resulting in a CV of the charge profile negative for highly charged cations and positive for highly 
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charged anions. For example in the case e+ > |e_|, e+p^ < |e_|pf* and q"* < 0. In this case {dip{z)/dz) > 0, 
J7+ > and the first term of eq. is positive. The second term in expression ^TE\\ is negative since for highly 
charged ions the coulombic part of the partial internal energy is more negative than for weakly charged ions and thus 
e+P+ + e_P_ < 0. However, this term must be sufficiently negative in order to partially cancel the first positive 
term in the expression of q'^^ and obtain an overall negative value for g*^'. For the case e+ < |e_|, g*^* will be positive. 
In this case, {dip{z)/dz) < 0, J+ < and the first term of g'^* is negative. The second term is positive and again it 
must cancel partially the first negative term in order to have g"^* > 0. We can see that the sign of g"^* for the charge 
asymmetrical case is the result of the partial compensation of the bulk term, which is the sum of the product of ionic 
charges with the partial pressure of ions, and the surface term connected with the electrostatic property J+. Both 
contributions have opposite signs but the bulk contribution must be dominant. 

The simultaneous consideration of charge and size asymmetry effects is more complex. Nevertheless, expressions 
(|18m9p for the CV of the charge profile are correct for this case also. It is clear that in some situations these two 
effects can enhance each others' influence and in other situations they can weaken or even cancel each other. 

In this paper, results previously obtained for symmetrical electrolytes @ are generalized for charge and/or size 
asymmetric systems. It is shown that the CV of the charge profile can be presented as the sum of three contributions. 
One of them is the normal component of the electrostatic Maxwell stress tensor. The second one is the electrostatic 
contribution from colons, described by the electrostatic property J7+ or depending on the sign of the surface. 
Similar as for the symmetrical case Q , J7+ or J7_ are defined by the integral from the product of the gradient of 
the electrical potential and the density distribution function of colons near the charged surface. The third term is 
the bulk contribution defined by the sum of the product of the charges of the colons and counterions multiplied 
by their corresponding partial pressures. Using the CTs for the charge and density profiles, the exact expressions 
for the CVs of the density profiles for colons and counterions are presented and discussed. We have also discussed 
the influence of the charge and size asymmetry on the CV of the charge profile for a noncharged surface. For 
noncharged surfaces, from the CT for the charge, the CV of the charge proflle now only has two contributions, 
the Maxwell tensor being zero. We show that in the case of size asymmetry both contributions have the same 
sign which corresponds to the sign of the CV of the charge proflle. In the case of charge asymmetry the two con- 
tributions have opposite signs. However the bulk contribution is dominant and imposes its sign to the CV of the charge. 

Finally, the generalization of the results we have obtained for the more than two components case faces the principal 
difficulty related with the problem of expressing p{z) using the charge profile q{z) and the colon density functions as 
in eq. pT|) . This problem can be solved only for the electrolyte with a single type of counterions. In this case, for 
example, for the positively charged surface 

q"*^ = ^e+aP+a{d+a/2) + e-p-{d-/2) 

a 

= /3e_— -^^e+„(e_-e+„)J+„+/3e_P_+^e+„P+„ (22) 

a a 

where 

J+a - / p+c.{z)^^^dz. (23) 

From eq. ([3]) and (|22|) we can obtain that 

^(e_ - e+a)p+a{d+/2) = /3^(e_ - e+a)P+Q + /3^(e_ - e+„)e+Q J+q. (24) 

However this equation is not sufficient to obtain the expressions for the CV of the ionic density profiles in a similar 
way as starting from expression p2p . 
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